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Abstract. We determine the average mass profile of galaxy clusters at two different redshifts and compare its evolution 
with cosmological model predictions. We use two samples of galaxy clusters spanning similar (evolutionary corrected) mass 
ranges at different redshifts. The sample of low-redshift (z ~ 0.0 — 0.1) clusters is extracted from the ESO Nearby Abell 
Cluster Survey (ENACS) catalog. The sample of high-redshift (z ~ 0.4 — 0.8) clusters is mostly made of clusters from the 
ESO Distant Cluster Survey (EDisCS). We determine the average mass-profiles for these two cluster samples by solving the 
Jeans equation for hydrostatic equilibrium, using galaxies as tracers. By using two cluster galaxy populations, characterized 
by the presence and, respectively, absence of emission-lines in their spectra ('ELGs' and 'nELGs' hereafter), we are able to 
partially break the mass-profile orbital-anisotropy degeneracy. 

We find that the mass-profiles of both the nearby and the distant clusters are reasonably well fitted by a Navarro, Frenk 
& White (NFW) model. The best-fit values of the NFW concentration parameter are as predicted by cosmological numerical 
simulations; cluster mass-density profiles become more concentrated with time. The evolution of the number-density profile 
of nELGs proceeds in the opposite sense, becoming less concentrated with time. 

In our analysis we also recover the orbital anisotropy of nELGs and ELGs . We find that in low-z clusters nELGs follow 
almost isotropic orbits and ELGs have more radially-elongated orbits. In high-z clusters both nELGs and ELGs follow radially- 
elongated orbits. 

We discuss these results in terms of the predicted secular mass growth of galaxy clusters and the transformation of ELGs 
into nELGs. 
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INTRODUCTION 

The study of the mass distribution within clusters and of its evolution with redshift can provide very useful constraints 
on the nature of dark matter (DM) HH^IH] and on the formation and evolution of galaxy clusters and their components 
. CDM cosmological simulations have shown that the mass-density profiles of cosmological halos follow a 
universal profile 0], 

Pnfw « (cr/r 20 o) _1 (l +cr/r 20 o) -2 , (1) 

parameterized by the concentration parameter c Q. This universal model, named 'NFW after the initials of its 
proposers, is characterized by a central density cusp. The NFW model has been shown to fit reasonably well the 
mass-density profile of nearby clusters [e. g . j | H I [Tol [Til [l2i [l3ll . although the precise form of the profile near the 
center has been a matter of debate [e.g. U4[l5l[l6ll . 

Cosmological simulations predict a mild dependence of c on the halo mass, M, and a mild evolution of the c—M 
relation with redshift, z re.g. ll7lll8ll . Observed c values of low-z galaxy systems appear to be somewhat higher than 
theoretical predictions, but the discrepancy mostly concerns low-mass galaxy systems lfl8l[l9ll . Massive galaxy clusters 
at low-z have the expected concentrations [e.g. [Toll20ll . This is true also for clusters at redshift z w 0.3 I2lll . 

Little is known about the mass-density profiles of galaxy clusters at still higher-z. Mass-density profile concentra- 
tions have so far been estimated for about a dozen clusters in total at z > 0.5 Il22ll23ll24ll25il26ll27il28ll29tl . all via 
gravitational lensing or deep X-ray observations. Here we report on the determination of the average mass-density 
profile of 19 clusters at 0.39 < z < 0.8 (z = 0.56), based on the projected phase-space distribution of cluster galaxies, 
used as tracers of the gravitational potential. 



The virial radius riaa is the radius within which the enclosed average mass density of a cluster is 200 times the critical density. 

The virial mass M200 is the mass enclosed within a sphere of radius r2oo. 

The circular velocity is defined from the previous two quantities as V2oo=(G A/200 / ^200 )'^ 2 - 



We adopt Hq = 70 km s 1 Mpc Q m = 0.3, £l\ — 0.7 throughout this paper. 

THE SAMPLES 

We use two samples of galaxy clusters, one at z = 0.56, the EDisCS sample ll30l[3ll[32i[33i[34ll . and another at z = 0.07, 
the ENACS sample I35LI36I1 . to investigate the evolution of the average cluster mass profile. 

We select the 15 EDisCS clusters with velocity dispersion (T/ as > 250 km s _1 , in order to have a more homogeneous 
data-set in terms of mass. We then add to this sample four clusters from the MORPHS ll3~7i l38tl . all with masses in 
the same range covered by the 15 EDisCS clusters. All 19 clusters have sufficiently wide spatial coverage (> 0.5r2oo) 
for the dynamical analysis, as well as homogeneous photometry (which is needed for the determination of the radial 
incompleteness, see below). 

Among the ENACS clusters, we use the 59 ENACS clusters studied in detail by l§9l [Tollioll . 

We identify cluster members by the procedure described in B41I1 . which has been validated on cluster-sized halos 
extracted from cosmological numerical simulations l4ll l42tl . We then determine cluster line-of-sight (los) velocity 
dispersions (T/ os by applying the robust biweight estimator to the velocity distributions of selected cluster members 
||43I1 . Cluster masses M200 are determined from the <7/ 01 -estimates using the 0/ o . s -M2oo relation of lliill . 

In summary, our data-set consist of 19 distant clusters from z = 0.39 to 0.80 (z = 0.56) with M200 masses from 0.7 
to 13.6 x 10 14 M Q (mean: 2.8 x 10 14 M Q ) and 59 nearby clusters from z = 0.03 to 0.10 (z = 0.07) with M 20 o masses 
from 0.4 to 20.5 x 10 14 M m (mean: 5.9 x 10 14 M Q ). When account is taken for the predicted evolution in mass between 
z = 0.56 and z = 0.07 11451 14611 the two samples are found to contain halos of similar (evolutionary-corrected) masses 
at two different cosmic epochs. 

Determination of the average cluster mass profile requires stacking together all the clusters in each of the two 
samples. Stacking is done by scaling the projected clustercentric galaxy distances, R, by cluster virial radii, r2oo, and 
the los galaxy velocities (in the cluster rest frame), v r / = (v — v)/(l +v/c), by cluster circular velocities, V2oo- In such 
a way we avoid mixing up the virialized regions of the more massive clusters with the unvirialized, external regions of 
the less massive ones. Since cluster mass profiles are expected to depend very mildly on cluster mass [e.g. 0,1471], the 
shape of the stacked cluster mass-profile is expected to be similar to the average shape of the individual cluster mass 
profiles. 

The center of each cluster is defined to be the position of its X-ray surface-brightness peak, when available, or the 
position of its brightest cluster galaxy otherwise. We only use galaxies in the radial range 0.05 < R/r2oo < 1 in the 
dynamical analysis. At R/rioo < 0.05 the analysis is unreliable because of the uncertainty in the position of the cluster 
center, and at R > r2oo the analysis is unreliable because dynamical relaxation is not guaranteed outside the virial 
region. In this radial range the high-z (low-z) stacked cluster contains 556 (respectively 2566) galaxies. 

We identify two populations of tracers: nELGs, the galaxies without emission lines in their spectra, and ELGs, 
galaxies with emission lines. More specifically, we classify ELGs the EDisCS galaxies with an [Oil] equivalent width 
> 3 A or with any other line in emission 13211 and the MORPHS galaxies with a spectral type different from 'k', 
'k+a', and 'a+k' IBal . We refer to 0511 for the ELG classification of ENACS galaxies. The fraction of ELGs among 
the galaxies selected for the dynamical analysis is 47% in the high-z sample, and only 13% in the low-z sample, a 
difference that reflects the evolution of the properties of cluster galaxies B37ll38ll32ll . 

The samples are not spectroscopically complete to a given magnitude. This is not a problem for the dynamical anal- 
ysis as far as the incompleteness does not depend on radius. The spectroscopic incompleteness is indeed independent 
on radius for the ENACS |36]. On the other hand, for the high-z cluster sample we need to weigh galaxies according 
to their radial positions when determining the galaxy number density profiles. The weighting method is described in 
Ml 

In constructing the galaxy number density profiles, we also need to assign different weights to the contributions of 
different clusters to the number counts in a given radial bin. This is because each cluster contributes galaxies to the 
stacked sample only out to a limiting radius, which is defined by the observational set-up. Hence, while at small radii 
all clusters contribute, at large radii we need to correct for those clusters that have not been sampled. The correction 
method is similar to the one described in fl48ll - see ll39ll and ifioll for applications to the ENACS. 




FIGURE 1. Left panel: The projected number density profiles, N(R), of the nELGs (red) and ELGs (blue) in the low-z stacked 
cluster. Right panel: same as left panel, but for the high-z stacked cluster. Solid lines represent best-fit models to the data, i.e. 
projected NFW models for nELGs and high-z ELGs, and the core model for low-z ELGs. 
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FIGURE 2. Left panel: The los velocity dispersion profiles Gi os (R), of the nELGs (red) and ELGs (blue) in the low-z stacked 
cluster. Right panel: same as left panel, but for the high-z stacked cluster. Solid lines represent best-fit models to the data, i.e. a 
c = 4.0 NFW mass model, a/rjrjrj = 3.6 (respectively a/rjoo = L2) OM velocity-anisotropy model for low-z nELGs (respectively, 
ELGs), and a.c — 3.2 NFW mass model, a/rioo = 0.01 ML velocity-anisotropy model for high-z nELGs and ELGs. 



THE METHOD 



The method we adopt for the dynamical analysis of our two stacked clusters is based on the standard spherically- 
symmetric Jeans analysis [490 In this analysis, the observables are the galaxy number density profile N(R) and the 
los velocity dispersion profile O[ 0S (R). These are displayed in Figures [T] and 12 respectively, separately for the two 
cluster samples and for the two cluster galaxy populations. 

N(R) is uniquely related to the 3-D galaxy number density profile v(r) via the Abel inversion equation. The 
other observable, &i os (R), is then u niquely determined if both the cluster mass profile, M(r), and the cluster velocity 
anisotropy profile, /3 (r), are known I50l 15111 . The velocity anisotropy profile /3 (r) is 

P( r ) = l ~^> (2) 



2 Spherical symmetry is imposed by the stacking procedure, which is irrespective of the individual cluster position angles 



where <v 2 >, <v 2 > are the mean squared tangential and radial velocity components, which reduce to of and a 2 , 
respectively, in the absence of bulk motions and net rotation (as we assume in the present analysis). 

We adopt parameterized model representations of M(r) and j3 (r) and determine the best-fit parameters of these 
models by comparing the observed &i os (R) profile with the predicted one, using the % 2 statistics and the uncertainties 
on the observed profile. In order to reduce the so-called "mass-anisotropy" degeneracy which plagues these kinds of 
analyses [see, e.g.,|52|,|2l],[IH] we adopt the method recently suggested by \ 5$\ . Namely, we consider two independent 
tracers of the same gravitational potential, nELGs and ELGs, and determine the best-fit parameters of M(r) and j3 (r) 
by a joint j 2 -analysis of the best-fits to the &i os (R) profiles of the two populations. Clearly M(r) must be the same for 
both tracers, but /3 (r) can in principle be different, so the degeneracy is only partially broken; however the constraints 
on the dynamics of the system are significantly stronger than when using a single tracer. 

Our choice of the M(r) and /3 (r) models is driven by the results of the analysis of cluster-sized halos extracted from 
cosmological numerical simulations. We adopt the NFW model, eq. (Till, parameterized by the concentration c, for the 
mass-density profile and also for the galaxy number-density profile but we allow different concentrations for the 
galaxy and the mass distributions. Only when the projected-NFW model does not provide an acceptable fit to N(R), 
we consider an alternative model, 

AToc [l + (R/R c f]- a . (3) 

We refer to this model|56] as the 'core' model, since it is characterized by a central constant density. 

We consider two models for the velocity-anisotropy profile j3 (r). One is the Mamon-Lokas ('ML' hereafter) model 
HI 

J8 = 0.5 r/(r + a), (4) 
and the other is the Osipkov-Merritt ('OM' hereafter) model ll58ll59tl 

P=r 2 /(r 2 + a 2 ), (5) 

Both the ML and the OM models depend on just one free parameter, the anisotropy radius, a, which marks the transition 
from the central region, where j3 w and the galaxy orbits are isotropic, to the external region, where j3 > and the 
galaxy orbits become increasingly radial. 



RESULTS 
The nearby cluster sample 

The nELG N(R) is best-fitted by a (projected) NFW profile with c = 2.4. A core model is required to fit the ELG 
N(R) which avoid the central cluster region, with best-fit parameter values Rc/r^ao = 1.28 and a = 3.2. The best-fit 
models are displayed in Figure \T\ (left panel). Abel-inversion of the N(R) best-fitting models provides the 3-D number 
density profiles v(r). 

The best-fit c parameter of the stacked cluster NFW mass-density profile is determined by a joint % 2 fit to the 
observed (7/ as -profiles of the nELGs and ELGs. The best-fit solution is obtained when the OM model is adopted for 
the velocity anisotropy profiles of the two galaxy populations. The best-fit value of the NFW concentration parameter 
is c = 4.0+1 3 (90% confidence levels, c.l. in the following), in agreement with [ 10]. The J 2 vs. c solution is displayed 
in Figure [3] (green curve). 

Using this solution for the cluster mass profile we obtain the best-fit j3 (r) OM-model parameters a/r2oo = 3.6+} 3 g 4 
and a / r2oo = 1 -2^q '\ for the nELG and ELG populations, respectively (90% c.l.). These j3 (r) profiles are shown in the 
left panel of Figure|4] [where we actually display <J r /<J t = (1 — j3)~ 1//2 ]. The best-fit j3 (r) solutions indicate that nELGs 
follow isotropic orbits within the cluster virial region, and that ELG orbits are isotropic near the center but become 
increasingly radial in the outer cluster regions. These results are in agreement with [40]. 

In Figure [2] (left panel) we display the observed los velocity dispersion profiles, and the best-fit obtained via the 
Jeans analysis. 



3 More precisely, we fit the projected NFW profile [ 55] to N(R). 
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FIGURE 3. The values of x 2 per degree of freedom obtained from the joint fit to the Oi os (R) profiles of nELGs and ELGs for 
different values of the c parameter of the NFW mass-profile model. The green and pink curves are for the low-z, and, respectively, the 
high-z stacked cluster. The two dashed lines indicate the 68% and 90% c.l. Symbols are theoretical predictions for the concentrations 
of the mass-density profiles of CDM halos with the same masses and redshifts as the clusters of our samples (squares: 1 17], dots: 
flUl : green: low-z, pink: high-z). The symbol error-bars are obtained by considering the 1-ct range in the distributions of cluster 
masses and average redshifts. 




FIGURE 4. Left panel: Best-fit velocity-anisotropy profile (7,-/(7, = (1 — /3)~ 1//2 for nELGs (red curve) and ELGs (blue curve) in 
the low-z stacked cluster. Dashed curves indicate 90% c.l. Right panel: same as left panel, but for the high-z stacked cluster. In this 
case the best-fit solutions are at the lower-limit of the interval considered in the x 2 minimization analysis, so only the lower c.l. to 
jS (r) are shown. 



The distant cluster sample 



Projected NFW models provide a good fit to the N(R) of nELGs and ELGs with best-fit parameters c = 7.5 and 
c = 2.7, respectively (see the right panel of FigureHJ. Abel-inversion of these best-fitting models then provide the 3-D 
number density profiles v(r) that we use in the dynamical analysis. 

The joint best-fit to the los (7/ ( „ -profiles of nELGs and ELGs is obtained when the ML j3(r) model is adopted. The 
mass profile best-fit concentration value is c = 3.2^ n (90% c.l.). The x 2 vs - c solution is displayed in Figure [3] (pink 
curve). Using this solution for the stacked cluster mass profile, we find that the best-fit value of the ML j3 (r) model 
parameter a, is the same for the nELGs and for the ELGs, a/rioo = 0.01, at the lower limit of the a-range considered 
in the % 2 minimization analysis. The solution is only poorly constrained, a < 0.9 r2oo and a < 10.0r2oo, for the nELGs 
and the ELGs, respectively. The best-fit velocity-anisotropy profiles are shown in Figure [4] (right panel). Both cluster 
galaxy populations are characterized by radially anisotropic orbits, although isotropy cannot be formally excluded for 
the ELGs, given the uncertainties. 

The observed o) OJ -profiles are shown in the right panel of Figure [2] together with the best-fit solutions from the 
Jeans analysis. 



DISCUSSION AND PERSPECTIVES 

Using cluster galaxies as tracers of the gravitational potential, we have solved the Jeans analysis for the dynamical 
equilibrium of a spherical system, and determined the average mass-density profiles of galaxy clusters at z — 0.1 
and at z — 0.6. These mass-density profiles are well described by NFW models with c values which decrease with 
increasing z. The best-fit c values are in agreement with the theoretical predictions of ACDM models. This can be 
seen in Figure [3] where we display our best-fit results together with the theoretically predicted mean c values. These 
are determined by using the masses and redshifts of the clusters in our two samples and applying the theoretical 
c = c{M,z) relations of lO] and IU8I1 . 

Our result appears to disagree with the conclusions of 1 18] who claim that the observed c values of clusters and 
groups are significantly above the theoretically expected ones for group- and cluster-sized cosmological halos. How- 
ever most of the claimed discrepancy is for low-mass galaxy systems, which are not covered by the present analysis. 
Preliminary dynamical analyses of the mass-density profiles of low-z groups 113], done with similar techniques as the 
one employed here, do indicate a discrepancy in the direction reported by 1181 . 

An additional result of our dynamical analysis is the determination of galaxy orbits in clusters at low-z and high-z. 
We find that these orbits become more isotropic with time. While low-z cluster nELGs have nearly isotropic orbits, 
high-z cluster nELGs move on radially elongated orbits, and so do both low-z and high-z ELGs (compare the left-hand 
and right-hand panels of Figure |4j. 

Orbital isotropization might result from the hierarchical accretion process of clusters [60] which undergo an initial, 
fast accretion phase, followed by a slower, smoother accretion pha se [46]. During the fast accretion phase clusters are 
subject to rapid variations of their gravitational potential JoH l62l l63ll . and these are capable of isotropizing galaxy 
orbits lo4l Iqj, l66l l67l l46Tl . The end of the fast accretion phase for cluster-sized halos occurs at z ~ 0.4 iPkill . hence it 
is not over yet for most of the clusters of our high-z sample. This can explain why the orbits of high-z cluster galaxies 
have not become isotropic yet. At lower-z, galaxies that have entered the cluster environment lately may experience 
only a slower orbital isotropization process, probably caused by ram-pressure [68]. Cluster ELGs are in this situation, 
while low-z cluster nELGs have developed isotropic orbits already because they were accreted much earlier than ELGs. 

The different relative fraction of nELGs and ELGs in low- and high-z clusters suggest that ELGs gradually transform 
into nELGs ll32ll . Since ELGs have a wider spatial distribution in clusters than nELGs, as ELGs join the nELG 
population, the global nELG spatial distribution become less concentrated. This is indeed observed in our data-sets 
(compare the red curves in the left and right panels of FigureQ]i. 

The results presented here are still preliminary, as their statistical significance is not very strong given the rather 
limited size of the data-set for high-z cluster galaxies. Substantial improvement in the understanding of the accretion 
and internal dynamics history of galaxy clusters requires much better statistics. In this sense, the proposed ES A mission 
"EUCLID'Q for the measurement of Dark Energy is extremely promising. As of this writing, the EUCLID mission 



4 http://sci.esa. int/science-e/ www/ area /index . cfm?f areaid=102 



proposal considers a 20,000 deg 2 spectroscopic survey with a near-IR slitless spectrometer down to a sensitivity level 
of 4 x 10~ 16 erg cm~ 2 s _1 for emission-lines, with a 1/3 sampling rate. Such a survey should return ~ 10,000 clusters, 
each with > 20 cluster members with measured redshifts, at 0.5 < z < 0.8. By stacking together clusters of similar 
mass it should be possible to constrain not only the evolution of the average cluster mass-density profile to an accuracy 
of <~ 3% in c, but also to measure the redshift evolution of the c = c(M) relation. With the same sample it will also be 
possible to determine the evolution of the ELG orbits in clusters, and to infer the history of cluster mass accretion. 
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